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Abstract — This paper investigates the performance of wireless 
systems that employ flnite-blocklength channel codes for trans- 
mission and operate under queueing constraints in the form 
of limitations on buffer overflow probabilities. A block fading 
model, in which fading stays constant in each coherence block 
and change independently between blocks, is considered. It is 
assumed that channel coding is performed over multiple coher- 
ence blocks. An approximate lower bound on the transmission 
rate is obtained from Feintein's Lemma. This lower bound is 
considered as the service rate and is incorporated into the 
effective capacity formulation, which characterizes the maximum 
constant arrival rate that can be supported under statistical 
queuing constraints. Performances of variable-rate and fixed- 
rate transmissions are studied. The optimum error probability 
for variable rate transmission and the optimum coding rate for 
fixed rate transmission are shown to be unique. Moreover, the 
tradeoff between the throughput and the number of blocks over 
which channel coding is performed is identified. 

I. Introduction 

An important class of wireless systems (e.g., systems that 
support streaming or interactive video applications) operate 
under quality of service (QoS) constraints in the form of 
limitations on queueing delays or buffer overflows. A measure 
of the throughput under such constraints is effective capacity 
[1], [2]. Effective capacity of wireless systems has recently 
been studied in various setting (see e.g., [l]-[4] and references 
therein). In most prior work, the service rates supported by 
the wireless channel are assumed to be equal to the instan- 
taneous channel capacity values and no decoding errors are 
considered. On the other hand, in practice, since finite block- 
length codewords are employed, communication is performed 
at rates less than the channel capacity with nonzero prob- 
ability of error. Recently, [5] has analyzed the performance 
of finite block-length codewords in the presence of statistical 
QoS constraints. However, in [5], coding is assumed to be 
performed over one coherence block in which the fading stays 
constant. In this paper, we consider a more general setting 
and assume that codewords are sent over multiple coherence 
blocks. Hence, each codeword experiences multiple fading 
realizations. Coding over multiple blocks generally improves 
the performance since codeword length can be increased and 
protection against severe fading can be provided as codewords 
see multiple channel states. However, coding over many blocks 
can also lead to long delays or buffer overflows. Therefore, it 
is of significant interest to analyze the throughput of channel 
coding over multiple coherence blocks in the presence of 
buffer constraints, and identify whether there exists an optimal 
number of blocks over which coding needs to be performed. 



II. System Model 

In this paper, we consider a block flat-fading channel, 
and assume that the fading coefficients stay constant for a 
coherence block of n symbols and change independently from 
one block to another. The discrete-time input and output 
relationship in the I th block is given by 



y l = hiXi +w { i = l,2, 



(1) 



where Xi and are the complex-valued channel input and 
output, respectively, in the i th symbol duration, hi is the 
channel fading coefficient in the Zth block, and W{ is the 
circularly symmetric complex Gaussian noise with zero mean 
and variance Nq, i.e., wt ~ CAf(0, Nq). We assume that the 
receiver has perfect channel side information (CSI) and hence 
perfectly knows the realizations of the fading coefficients {hi}. 
On the other hand, we consider both cases of perfect and no 
CSI at the transmitter. 

The channel input is assumed to be subject to E{|xi| 2 } < 
£ s . It is well-known that when the receiver has perfect CSI, 
the capacity achieving input for the above fading Gaussian 
channel is Gaussian distributed. Hence, we assume that xi ~ 
C7V(0, £ s ). Since the input and noise are Gaussian distributed, 
the output is also conditionally Gaussian, i.e., yi\hi ~ 
CN{Q,£ s \hi\ 2 + N ). Moreover, yi\xi,hi ~ CN{h lXl , N ). 
We further assume that the input is independent and iden- 
tically distributed (i.i.d.) i.e., p x ™ = Y\2=iPxi( x i)> which 
implies p y n\ x n thl = n"=i Py t \x it h, (Vi\%i, hi), and p y n\ hl = 

U7=iPy t \h l (yi\ h i)- 

III. Preliminaries 

A. Effective Capacity 

In [2], effective capacity is defined as the maximum constant 
arrival rat^that a given service process can support in order to 
guarantee a statistical QoS requirement specified by the QoS 
exponent 8. If we define Q as the stationary queue length, then 
9 is the decay rate of the tail distribution of the queue length 



Q: 



lim 

q— >oo 



logP(Q>g) 



q 



(2) 



Therefore, for large <7 max , we have the following approxima- 
tion for the buffer violation probability: P(Q > q max ) f» 
e -e?max Hence, while larger 9 corresponds to more strict 
QoS constraints, smaller 9 implies looser QoS guarantees. 

For time- varying arrival rates, effective capacity specifies the effective 
bandwidth of the arrival process that can be supported by the channel. 



Similarly, if D denotes the steady-state delay experienced in 



the buffer, then P(D > d max ) ~ e 



-8Sd„ 



for large d„ 



where S is determined by the arrival and service processes 
[3]. The effective capacity is given by 

Re{0) = - lim — log e E{e- es W} bits/s, (3) 

where the expectation is with respect to S[t] = X)!=i r *M' 
which is the time-accumulated service process. {r s [i],i = 
1,2,...} denotes the discrete-time stationary and ergodic 
stochastic service process. 

B. Mutual Information Density and Channel Coding Rate 

As detailed above, effective capacity is determined by 
specifying the service rate or equivalently the instantaneous 
transmission rate. We assume that the transmitter performs 
channel coding over m coherence blocks where m = 1,2,... 
Therefore, it sends codewords of length nm and each code- 
word experiences m independent channel conditions. An upper 
bound on the maximum decoding error probabilities of random 
codes of length nm is given by Feinstein's Lemma [6], [7]: 



, nm 



+ P(x nm i S nm ) + e 



(4) 



where 7 > is an arbitrary constant, S nm = 
{^Er=TiE{M 2 }<£} is the constraint set, 
i(x nm ; y nm \h™) is the mutual information density conditioned 
on the fading coefficients (hi, hi, . . . , h m ) seen in m 
coherence blocks. The conditional mutual information density 
is defined as 

„/,.nmLnra h m \ 

i(x ;y \h 1 )=log 2 — 7 ^———. (5) 



p(y 

Next, we obtain an expression for the mutual information 
density of the considered channel and input models (i.e., fading 
Gaussian channel with Gaussian input), and derive, under 
certain assumptions, an approximate lower bound on the rates 
attained by coding over m coherence blocks. 

For the system model introduced in Section HU , we have 

1 . / nm nm i 7 m \ 

i(x ;y \hi ) 



^ m In 

— i(xi;yi\hi 



1 = 1 i = (Z-l)n + l 
m In 



= — E E lo ? 

nm z — ' z — ' 



fy^htiVi, hi) 



1=1 i=(I-l)«+l 
m In 



;=i i=(;_i)n+l 
|2 



£s\h 



Vi\ l°g 2 e — /ijacil log 2 e 



\hi\ 2 £ s +N 



No 



1=1 



£4M 
iVo 



; = 1 i=(!-l)n+l 



|/u| 2 £ s + iVo No 



Denoting SNR = jf- and extending the results in [6] and 
[7], we can immediately show that i(x nm ; y nm ) / (nm) has 
the same distribution as the random variable [8] 



1 JUL l 

lE 1 °g 2 ( 1 + SNR i /i 'i 2 ) + ^ 

in — ' nm 



E 



SNR|fe;| : 



1=1 



nm ^ 1/ 1 + SNR|/i ; | 2 ^ 
i=i V 1 " i=i 



?/'/, 



where iu^'s are i.i.d. Laplace random variables, each with 
zero mean and variance 2. The sum of nm i.i.d. Laplace 
random variables has a Bessel-K distribution [6] and generally 
is difficult to deal with directly. On the other hand, for 
large enough values of the blocklength nm, it can be well 
approximated by a Gaussian random variable [7]. Therefore, 
the mutual information density achieved with the codewords 
of length nm spreading over m coherence blocks can be 
approximated as 

— i(x nm ;y nm ) 



nm 
~CJV 



— > log 2 (l + SNRZ/), > 



(6) 



where we have defined zi = \hi\ 2 . With this approximation, 
the first probability expression on the right-hand side of © 
can be written in terms of the Gaussian Q-function: 



P I —i(x 
nm 



,nn -y nm \hT) <R + i 



= Q 



^Etll0g 2 (l + SNRz i )-i?-7 



log| e 



El 



2SNR.S 



(7) 



A=l nm(l+SNRz,) 

By noting that the Q-function is invertible, we can rewrite the 
upper bound in l|4} as a lower bound on the instantaneous rate 
achieved by coding over m coherence blocks: 



R >— V log 2 (l + SNRz ; ) 

m z — ' 



i=i 



\ 



E 



2SNRZ; 



71771 (1 + SNRZ; ) 



x Q-\e - P(x nm i S nm ) - e- nm ~<) - 7 



(8) 



for any 7 > 0. Although the above lower bound can also be 
used in the subsequent analysis, we opt to further simplify it to 
make the analysis more tractable analytically. For sufficiently 
large values of nm, the terms P(x nm ^ S nm ) and e - "™ 7 
become very small and can be neglectecQ Moreover, since 
the lower bound holds for any 7 > 0, we can see that an 
approximate lower bound for the transmission rate is 



R > Rl 



m 

-E l0 S2(i 

m * — » 



SNRZH 



1=1 



\ 



2SNRZ/ 



nm(l 



SNRZ;) 



Q-'ie) (9) 



2 As nm increases without bound, it can be easily seen that P(x nm £ 
Snm) approaches zero by noting the fact that the codewords are generated 
according to p x nm = f]™^))^ (x;) where p Xi is the Gaussian distribution 
with zero mean and variance £ and by applying the law of large numbers 
which tells us that the sample variance approaches the statistical variance £. 
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where the notation e is used to emphasize that this is a 
lower bound for rates achieved with decoding error probability 
e. Henceforth, the analysis is based on 

IV. Effective Throughput with Channel Coding 
over Multiple Coherence Blocks 

The rate lower bound in (|9jl gives a characterization of the 
tradeoffs and interactions between the instantaneous transmis- 
sion rate, decoding error probability and the fading coefficients 
when channel coding is performed over multiple coherence 
blocks using finite blocklength codes. In particular, we note 
that Ri e is achieved with probability 1 — e. With probability e, 
decoding error occurs. We assume that the receiver reliably de- 
tect the errors, and apply a simple ARQ mechanism and sends 
a negative acknowledgement requesting the retransmission of 
the message in case of an erroneous reception. Therefore, the 
data rate is effectively zero when error occurs. Under this 
assumption, the service rate (in bits per nm symbols) is 

fO, with probability e 

s y nmRi if , with probability 1 — e 

Similarly as in [5], we obtain the following result on the 
effective rate by inserting the above service rate formulation 
into the definition in Q and noting that the service rate varies 
independently for one sequence of to blocks to another due 
to the block fading assumption. Since it characterizes the 
throughput achieved by transmitting at rates possibly below 
the channel capacity using finite blocklength codes, we refer to 
this throughput measure as effective rate rather than effective 
capacity in the remainder of the paper. 

Theorem 1: The effective rate (in bits per channel use) at a 
given SNR, error probability e, codeword length nm, and QoS 
exponent 9 is 

Re(9) = --Llog e E B {e + (1 - e)e- 6nmR ' <} (11) 

where Ri. e is given in (|9j, and z = (z±, . . . , z m ) is the vector 
composed of the channel states experienced in m blocks. 

The effective rate in (II Q provides a lower bound on the 
throughput as a function of SNR, decoding error probability 
e, fading coefficients, coherence blocklength n, the number 
of blocks, to, over which coding is performed, and the QoS 
exponent 9. The following result shows that given the other 
parameters, the effective rate is maximized at a unique decod- 
ing error probability. Note that using very strong codes and 
having small error probabilities in the transmission necessitates 
small transmission rates leading to small throughput. On the 
other hand, if higher transmission rates with relatively weak 
channel coding are preferred, then communication reliability 
degrades and more retransmissions are required again lowering 
the throughput. 

Theorem 2: Given the values of m > 0, n > 0, 9 > and 
SNR > 0, the function 

*(e) =E{e + (l-e)e- 0nmR '''} (12) 

is strictly convex in e and hence the optimal e > that mini- 
mize ty(e), or equivalently maximizes the effective throughput, 
is unique. 




Fig. 1. The effective rate as a function of e. n = 200. 



Proof: See Appendix A. 

In Fig. ^ we pl°t the effective rate as a function of e in 
the Rayleigh fading channel with E{z} = 1. In the figure, we 
assume SNR = dB, 9 = 0.01, and n = 200. We provided 
curves for different values of to. We can see that the effective 
throughput is indeed maximized at a unique e. 

Another interesting tradeoff is the one between the through- 
put and the number of block over which coding takes place. 
Note that coding over multiple blocks is generally beneficial 
in terms of increasing the transmission rate on the average 
because transmitted codewords experience multiple channel 
fading realizations and may not get exceedingly affected by 
severe fading in one block. On the other hand, coding (with 
possible retransmissions in cases of decoding errors) over 
multiple blocks may lead to unacceptable delays in systems 
operating under buffer constraints captured by the QoS ex- 
ponent 9 in this paper. Therefore, coding over many blocks 
may lead to decreased throughput eventually. In Fig. [2] we 
plot the effective throughput as a function of to for different 9 
values with fixed e = 0.01. We assume n = 50. In the figure, 
we observe that the optimal m that maximizes the effective 
rate under a given e varies with 9. When 9 — and therefore 
there are no buffer constraints, effective rate increases with 
increasing to. Coding over ever increasing number of blocks 
improves the performance. Indeed, as to — > oo, effective rate 
approaches the ergodic capacity in the case of 9 = 0. However, 
we see a strikingly different behavior in the presence of QoS 
limitations. We note that for 9 > 0, effective rate is maximized 
at a finite value of to. Moreover, the optimal value of to 
diminishes as 9 increases. Therefore, coding over fewer blocks 
should be preferred under stringent buffer limitations. 

In Fig. 13 we plot the optimal effective rate (optimized over 
the decoding error probability, e) as a function of 9 for given 
to values. We set n = 50. We find that for small 9 values, 
having to = 10 achieves the highest effective rate, while as 9 
increases, having to = 10 starts providing the lowest effective 
throughput, due to similar reasons as outlined above. 

A. Fixed Rate Transmissions 

Heretofore, we have implicitly assumed that the transmit- 
ter has perfect CSI and considered the scenario in which 
the transmitter employs variable-rate transmissions with rates 
characterized by Ri £ given in l|9}. Note that in order to 
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Fig. 2. The effective rate as a function of m. n = 50. e = 0.01. 



Fig. 4. The effective rate as a function of R. n = 200. 




SNR > are fixed, then the function 



Fig. 3. The optimal effective rate vs. 8. n = 50. 



transmit at the rate i?j j£ , the transmitter needs to know the 
fading coefficients {zz}™^ A more practical scenario is the 
one in which the transmitter does not know the channel states 
and send the information at a fixed rate of R. Note that 
in this case, the decoding error probability varies with the 
fading coefficients in each set of m blocks in contrast to 
being constant in the previous analysis. The codeword error 
probability for a given channel state z is 



e{z,R)=Q 



m 1^1=1 



2SNRz 



(13) 



nm(l+SNR2() 1°§2 e , 

obtained by using (|9}- The effective rate is then expressed as 
Re(0,R) 

1 



mm 
1 

9nm 



log e E z {e(z, J R) + (l-e(z,J?)> 



-9nmR 



} 



■logMQ 



1-Q 



/ 1 v^m 2SNR2, i 

m 1^1=1 nm(l+SNRz,) iU &2 e 

^E;"llQg 2 (l + SNRzQ-^ 

AT v^ m 

m Z^i=l 



— QnmR 



nm(l+SNR 2i ) 10 §2 e , 



(14) 



We have the following result. 

Theorem 3: Assume that the values of n, m, 9 > and 




• 2SNRz, i 

= 1 nm(l+SNR Z! ) iU &2 1 



^E£ilog 2 (l + SNR«l)-^ 



-QnmR 



m Z^Z — 1 



2SNRz ; 
nm(l+SNR2( 



log 2 e 



is minimized at a unique i? and hence the optimal R that 
minimizes Q(R) or equivalently maximizes the effective rate 
in is unique. 
Proof: See Appendix B. 

In Fig. |4] we plot the effective rate as a function of the 
fixed transmission rate R for different m values. We assume 
that 8 = 0.01, n = 200. It is noted that there is a unique 
R that maximizes the effective throughput. While we see in 
the figure that the optimal effective throughput increases as 
m increases from 1 to 10, further numerical analysis showed 
results similar to those discussed in the variable-rate case, i.e., 
effective rate starts decreasing when m is increased beyond a 
threshold value. 

V. Conclusion 

We have analyzed the performance of channel coding over 
multiple blocks with possible decoding errors in the presence 
of queueing constraints. We have characterized the effective 
throughput. We have discussed two different transmission 
strategies. For the case when the transmission rate is varied and 
the error probability is kept fixed over different codewords, we 
have shown that the optimal error probability that maximizes 
the effective throughput is unique. Similarly, when the trans- 
mission rate is kept fixed, we have proved that the optimal 
rate that maximizes the effective throughput is also unique. 
Through numerical analysis, we have quantified the tradeoff 
between the throughput and the number of blocks over which 
coding is performed. 

Appendix 

A. Proof to Theorem 

We first prove the following. 
Proposition 1: The function 

/(e) = (1 - e )e- 6nmRl ^ 



(15) 



4 



is strictly convex in e. 



Proof: Denote 

-OnrnRi, e = aQ^ 1 (e) + b 
where we, from l|9}, define 



(16) 



— — — log 2 e, b = -On ^ log 2 (l + SNR*,). 



(17) 



Note that a > since SNR > 0, > and > 0. Then we 
can rewrite (I15> as 

/(e) = (l-e)e aQ ~ 1(e)+b . (18) 

The first and second derivative of /(e) with respect to e are 

/(e) = (aQ- 1 (e)(l - e) - l) e^M+i (19) 

/(e) = (a(l - e) (<r x (e)) 2 - 2Q- 1 (e) + (1 - e)Q- 1 (e) 

x ae aQ_1(e)+b (20) 

where Q _1 (e) and Q _1 (e) denote the first and second deriva- 
tives of Q _1 (e) with respect to e. Note that for an invertible 
and differentiable function g, we have g(g~ 1 (x)) = x. Taking 
the derivative of both sides, we have 

yig-'ix^g-^x) = 1 g-^x) = 1 (21) 

9(9 1 (x)) 

where g(g~ 1 (x)) denotes the first derivative of g evaluated at 
g^ 1 (x), and g^ 1 (x) is the derivative of g^ 1 with respect to 
x. Noting that 



Q(x) 



e 2 dt, Q(x) 



e 2 



we can derive the following 



2ire 



(22) 



(23) 



Note that Q~ l {e) < for < e < 1. Differentiating Q -1 (e) 
with respect to e, we have 

Q- 1 (e) = 2 7 rQ- 1 (e) e -^ 1 ( £ )) 2 . 



(24) 



Next, we consider the following two cases: 



1) e < |: We have Q^ 1 (e) > for this case and hence 
Q _1 (e) > 0. Together with the fact that Q" 1 ^) < 0, we can 
immediately see that /(e) > for e < |. 



2) e > ±: We have Q~ l {e) < for this case. Substituting 
d23l and (1241 to (I20> and denoting a; = Q _1 (e), the expression 



inside the parentheses of (I20> can be written as 

o(l - e) fQ-^e)) 2 - 2Q" 1 (e) + (1 - e)Q~\e) 



= o(l - e)27re" (Q 1(e))2 + 2\/27re 

+ (l-e)27rQ- 1 (e)e- crl(e))2 



(25) 



(26) 



= o(l - Q(x))27re"— + 2V27r e — + (1 - Q(x))2wxe~ x 

(27) 



(2vr(l - Q{x))(x + a)e^ + 2\/2n S j 
(2n(l - Q(x))xe s r + 2V2tt) 



> 

> eT" I 2?r 



-( - .;v 2 + 2V27T 

/ 2^> 



e 2 



(28) 
(29) 
(30) 
(31) 



/2vr(-x) 

= (-%/2tt + 2%/2tt 

where ( I29l l follows from the facts that a > and hence aj+a > 
x. d30l is obtained from the following upper bound 

1 



1 - Q{x) = Q(-x) > 



-e 2 for x < 



(32) 



/2tt(-^) 

and the fact that x < for this case, and hence multiplication 
of x(l — Q(x)) can be lowerbounded. Therefore, /(e) > for 



e > 



Also note that e = A means Q 1 (e) = 0, so we have 



a(l - e) (Q-^e)) - 2Q-\e) + (1 - e)Q~ 1 (e) 

= o(l - e)2^e- W( £ )) 2 + 2V2^e- i ^^ 

+ (l-e)2^Q- 1 (e)e- <3 "( £ » 2 
= a?r + 2v / 2?r > 



(33) 

(34) 
(35) 



and as a result /(e) > 0. 

From the above discussion, we can find that /(e) > for 
all e € [0, 1]. /(e) is strictly convex in e. □ 

Now, let ijj(e) = e + (1 - e ) e - 6nmR ^ = e + /(e). 
ip(e) — /(e) > 0. Since the nonnegative weighted sum of 
strictly convex functions is a convex function [9], we can 
conclude that ^(e) is convex. 



B. Proof to Theorem 

First, for any given channel state pair z = (zi, Z2, ■ ■ • , z m ), 
we define 



= - Vlog^l + SNR^), 



(36) 



i=i 



5 = 



\ 



1 ^ 2SNRZ; 

— > t, r lo §2 e (37) 

m nm(l + SNRz;) 



and note that // > 0, (5 > 0. We can find that $(0) = 
1, $(oo) = 1, and $(i?) < 1 for all R £ (0,oo). Note that 



Q(x) = 



1 



(38) 



5 



The first and second derivatives of &(R) in R are given by 
-0nm(l-E{Q(^^]\)e- 9nmH <»> 



2 (f-R) 2 i / f fu — R 

x | 1M : e 24a ^ + fa™ 1 - E ' ' ' ' '" 

2-kS 



(40) 



Now we need the following result. 

Proposition 2: <l>(i?) = has only one solution. 
Proof: Obviously, $(0) > 0. Letting <J>(i?) = and perform- 
ing a simple computation, we have 



E 



/2ttS 



E 



9nm 



1 — E 



6nm- 



V 



E 



9nmR ' 



(41) 

First, we can show that the left-hand side (LHS) of J41i is an 
increasing function in R. Let 

1 



= E 



2tt<5 



(42) 



i 

. c e j s log-concave function [91. Since the nonneg- 

V27TO 

ative weighted sum of concave functions is concave, g(R) is 
log-concave function [9]. And hence — log e <7(i?) is convex 
function. Note that 

LHS = ^(-log e 9(R)) (43) 

thus the derivative of LHS of ( 14 1> is greater than 0, and as a 
result it is increasing in R. 

Next, we can prove that the right-hand side (RHS) of (1411 is 
a decreasing function in R. Note that ^_ Q nmH is decreasing 
in R. Let 

* v ( u — R \ 

^ji R "■ — 



1- 



Q 



= E- 



1 

/2jr 



hence J41i has only one solution. □ 
Denote the unique solution to <J>(i?) = as R' . We know 
that <3?(i?) > for all R < R', or <f>(i?) is increasing 
equivalently, and < for all R > R', or &(R) 

is decreasing equivalently. Note here that J °° &(R)dR = 
$(oo) - $(0) = 0, 6(0) = -6nm{l - E{Q(§)}) < 0, 
so <l>(i?') > 0. Otherwise, $(i?) is decreasing for i? > i?', 
and hence < 0, /°° §{R)dR < 0. Contradiction. Also 

note that $(oo) = 0, so $(i?) > for i? > i?'. Thus, there 
is only solution to <P(i?) = in the range R g (0, i?'), and 
is minimized at this value. 
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e-Trftj =g(u(R)), 
(44) 



We 



where g(x) = E j J*^ ~^^ e * 2 dt\ and u(x) 
know that g(x) is a log-concave function [9], and from [9, 
Eq. 3.10], we can see that \og e g is concave function, and u 
is concave and nonincreasing in R, and hence log e g(u(R)) 

E 

is concave function in R directly.So — 



q»— «)- 

94^ 



_E{Q(^)| 
l-E{Q(ii^)} 

increasing function, i.e., — — ^-r^ \ ls a decreasing 



function in R. Thus, the RHS of <!41i is decreasing in R, and 
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